arXiv: 1502.07711v3 [hep-th] 23 Apr 2015 


Scattering equations, generating functions and all 
massless five point tree amplitudes 


Chrysostomos 


Kalousios,“ 0 


“ ICTP South American Institute for Fundamental Research 
Instituto de Fisica Teorica, UNESP-Universidade Estadual Paulista 
R. Dr. Bento T. Ferraz 271 - Bl. II, 01 If 0-070, Sdo Paulo, SP, Brasil 


Abstract 

We argue that one does not need to know the explicit solutions of the scattering equations 
in order to evaluate a given amplitude. We consider the most general quantity consis¬ 
tent with SL{2, C) invariance that can appear in an amplitude that admits a scattering 
equation description. This quantity depends on all cross ratios that can be formed from 
n points and we evaluate it for the first non-trivial case of n = 5. The combinatorial 
nature of the problem is captured through the construction of an appropriate generating 
function that depends on five variables. 
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1 Introduction 


In [i| it was argued that the tree level S-matrix of massless theories can be captured by 
the so-called scattering equations (to be dehned later) that connect the space of kinematic 
invariants of n particles in arbitrary spacetime dimensions to the positions of n points on 
a sphere. Prior to [H the scattering equations had appeared in the literature in different 
contexts in 0-[i3. 

After the initial conception and application of the scattering equations to Yang-Mills 
and gravity llj], an increasing number of theories, whose tree level amplitudes can be 
expressed in terms of the scattering equations, has been found [iMi some of them 
generalizing to massive cases, with the promise that this is not the complete list. It is still 
not known what kind of theories admit a representation in terms of scattering equations. 

The formalism was proven for Yang-Mills in 1^, where the authors showed that 


it reproduces the BCFW [l9j recursion relations. 


equations that greatly facilitates computations was also presented in 


A polynomial form of the scattering 

Extensions at 


loop level include [2l|-|23l| and connection to twistor-string-like models can be found in 


24-26 


One of the questions is how to use the formalism in order to get explicit answers for 
the amplitudes. Some attempts have appeared in the pas t and involved special solutions 
of the equations associated to particular polymonials Il2l . l20l . I27H29I , that in some of the 


cases allowed for the explicit construction of the amplitude '|l2l. l27j|. Besides all the efforts 
there is no known general solution of the scattering equations. 

Fortunately, one does not need to know the explicit solutions of the equations in order 
to evaluate the amplitude. The amplitude is always given as sums of all possible solutions 
of the scattering equations, which are polynomial in nature, and one can then use the well 
known in mathematics formulas of Vieta, that associate the sums of roots of polynomials 
to the coefficients of these polynomials. This can be the hrst step, but still the answer is 
complicated to evaluate and write it down. 

In this work we attempt to organize the expressions of tree level amplitudes. We hrst 
identify a fundamental quantity, whose general form depends on products of cross ratios. 
We claim that all amplitudes can be written as linear combinations of this quantity. 
In the case of n = 5, we explicitly evaluate this fundamental quantity by constructing a 
generating function that captures the combinatorics of the problem. We then give specihc 
examples for the case of Yang-Mills. 


2 Scattering equation formalism 


The fundamental ingredient that allows the formulation of S-matrices in arbitrary dimen¬ 
sions is the scattering equations and are dehned as 


/. = Y 

h^a 


^ab 

^ab 


( 1 ) 


where ka is the momentum of the particle. In the above we have used the short 
notation kab = ka-kb and aab = era — Cb- Not all of the n equations in ([1]) are independent. 


1 
































but instead they satisfy three constraints 




This is due to the SL{2,C) invariance of ([I]), which is a direct consequence of total 
momentum conservation and the on-shell condition of the external particles. This allows 
us to £x three of the cXjS to arbitrary values. The number of solutions of ([T]) is known [l| 
to be (u — 3)! and in general they can be complex. 

The S-matrices of the theories that admit a scattering equation description have the 
general form 

/ d’^CT -I—r 

vol SL(2 C) n 

where In{k,e,a) depends on the theory and carries information about the external par¬ 
ticles, namely their momentum k and polarization vectors e. Invariance of the integrand 
in (|3]) under S'L(2,C) transformations restricts the form of /„(fc,e, a). Finally, the delta 
functions appearing in ([3]) completely localize all integrals. 

We will now present how the scattering equations are used in the case of Yang-Mills 
amplitudes, since later we will make use of the formulas in our examples. From [lit] we 
know that after performing the integration ([3]) the tree level n-gluon partial amplitude 
An of Yang-Mills in arbitrary dimensions can be expressed through 


A 


n 


E 

roots 


1 PhT(fc,6,cr) 

cri2a23 ■■■(Tni det'$ 


(4) 


where the sum runs over all solutions of ([T]). The 2n x 2n antisymmetric matrix T is 
given by 

C B 


T = 


where the n x n matrices A, B, C are given by 

^a'^b e 


ka'k}) 

^ab — ^a^b-) 

(^ab 


Bnh — 






^ ^a'kb r- 

^ab Oa^b 

CTab 


Sab 


en-kr 


Cy^a 


(Jn 


The matrix $ is dehned as 


dfa ka-kb -^ka-kc 

^ab = — - = - —Oajtb — Oab 2_^ -^ 


dab 


a. 


ab 


c^a 


(Jn 


and the primes in (jl]) denote 


Pf'T = 2 ^ det'd) = 


det(<F;^J) 


a., 






rp J 


(5) 


( 6 ) 


(7) 


( 8 ) 


th 


The matrix in ([8]) is derived from the matrix T after the removal of the and j 
row and the and column, with 1 < i < j < n. Finally, the matrix is derived 
from the matrix <F after removing the {q j, k} rows and the {p, g, r} columns. 
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3 An algorithm to evaluate the amplitudes 


The first question we would like to address is whether the scattering equations can be 
explicitly solved and whether such a solution is useful for calculating the amplitudes. As 
it was shown in 20j the scattering equations admit a polynomial form. According to 


one can use the elimination theory and completely decouple the scattering equations. 
Then one ends up with a one variable (n — 3)! degree polynomial equation, p((Jj) = 0, 
for one of the (n — 3) variables Uj, whereas the rest (n — 4) variables can be uniquely 
determined from the solution of the aforementioned polynomial. The coefficients of p{(Ji) 
depend only on products and sums of the kinematic invariants ki-kj and in general can be 
very long to write them explicitly. For the simplest non trivial case, n = 5, the solution 
of p{<Ji) = 0 already occupies several lines, whereas for the next case, n = 6, as explicit 
calculations for special kinematics show, the six solutions become long in an uninspiring 
way and extend into several pages. 

In general it is not known whether one can explicitly solve any higher case. Fortunately, 
one does not have to do so in order to evaluate the amplitude. The explanation is simple. 
One can in principle obtain all amplitudes of any known theory that admits a scattering 
equation description using the following procedure, which does not require the explicit 
solution of the scattering equations. The idea is the following. One uses the results 
of in order to write the amplitude in terms of only one variable. The answer for 
the amplitude will then be a ratio of two polynomials of the same variable of degree 
much higher than (n — 3)!. Then one can iteratively use the scattering equation of the 
remaining variable and bring the amplitude to the form of a ratio of two polynomials of 
degree (n — 3)! — 1. Then one can use the well known Vieta formulas that associate the 
sum of roots of a polynomial to its coefficients and obtain the amplitude as a rational 
function of the kinematic invariants. 

Let us illustrate the algorithm with the help of the following toy model. We consider 
the toy scattering equation to be 


— ax + b = 0 


and the toy amplitude to be 

4 

roots 

Iterative use of ([9]) in flTOD gives 


^toy — 

roots 


(a^ — 2ab)x + (6^ — a% — c) 
(a^ — h)x — {ab + d) 


(9) 

( 10 ) 


( 11 ) 


Let the two solutions of ([9]) be denoted by ri and r 2 - Substitution to flTT]) yields 

(a^ — 2ab)ri + (6^ — a‘^b — c) {a^ — 2 ab)r 2 + (6^ — a‘^b — c) 
(a^ — b)ri — {ab + d) (a^ — b)r 2 — {ab + d) 

_ Cl + C 2 (ri + r2) + C3rir2 _ Ci + C2a + C 36 
C 4 + C5(ri + r2) + C6rir2 Ci + c^a + Cob' 


( 12 ) 


where the constants Cj are simple functions of a, 6, c, d that can be easily computed. In 
the above we have made use of the Vieta formula that states that the sum of roots of 
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([9]) is ri + r 2 = a and the product rir 2 = h. Our toy model can be easily applied to 
the most general case. One has a polynomial scattering equation of degree higher than 
two, therefore the hnal expression of the amplitude contains not only the sum and the 
product of roots of the scattering equation, but all the elementary symmetric polynomials 
of the roots. Hence, we conclude that one does not need to solve the scattering equations, 
whereas at the same time we have shown that the amplitude is a rational function of the 
kinematic invariants ki-kj as expected. 

Although the above algorithm does not require the explicit solution of the scattering 
equations, it becomes quickly complicated. The authors of have stopped the demon¬ 
stration of their construction at n = 6 . Although in principle it is possible to extend 
the analysis to higher cases, it becomes difficult to continue beyond n = 6 or n = 7 and 
one should perhaps rely on other ideas in order to explicitly obtain, organize or write 
down the amplitude. One such idea is the expression of the amplitudes with the help of 
a generating function, to which we now turn. 


4 Calculation of the generating function 

All n = 5 amplitudes can be decomposed as sums of the following fundamental quantity 


Pa = 


roots 


1 1 
det'<h ai2^+"la23^+"2o-342+«3o-452+a4o-^52+a5cr^3a6o-i4«7cr24«8o-25"9cr35“l0 ’ 


(13) 

with momentum and helicity dependent coefficients. In the above, the ttjS are assumed 
to be integers. Such a decomposition might not be immediately obvious and we explain 
how this can be done later in an example. We demand SL{2, C) invariance that will £x 
hve of the ajS to the values 


OlQ — —Oi — 0.2 T O 4 , 
as = ~Q^2 ~ 0^3 -f- as, 
aio = Tai — a3 — a4. 


Oj — “l~a2 — a4 — as, 
Og = —Oi -f- a3 — as. 


(14) 


In the dehnition of P in flT^ one can shift any of the a^s by an integer value without 
loss of generality. Since there is no canonical way to express that, we have chosen P to 
correspond to the color ordered 0^ amplitude jl^ when we set a* = 0, z = 1,..., 5. Other 
starting points can also be considered and lead to hnal expressions of the same or higher 
complexity. 

Substitution of ffTTj) in (IT^ yields 


roots 


1 

det'd) 


n 


i=l 




o^i 


( 15 ) 


The number of cross ratios appearing in flT^ coincides with the number of independent 
cross ratios in d-dimensions, namely n{n — 3)/2. Since our problem is one dimensional the 
number of independent cross ratios in our case is n — 3, which means that the conformal 
ratios in (lT5ll are dependent, in general in a complicated way. We have traded away this 
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complication by considering a larger set of cross ratios that has the advantage that all 
remaining cross ratios can be simply expressed as products of quantities of that set. 

The scattering equations for the n = 5 case are quadratic in nature and admit two 
solutions that contain one square root. After we £x the SL{2, C) invariance and substitute 
the solution of the scattering equations in P we get an expression of the following form 

P = ^{bo + coy/r) + hbo - co^/r) - Ci^/rY\ (16) 


where the foo, cq, Q, r appearing in the above expression are rational functions of the 
kinematic invariants ki-kj and in general can be complicated. bo,co depend on how we 
choose to parametrize the ctj = 0 case. For integer a^s the quantity P is a rational function 
of the kinematic invariants. Here and in the rest of this work the index i always takes the 
values 1,..., 5 and never zero. We consider cyclicity of indices and identify i + 5 ~ i 
There is a nice way in mathematics to express flT6|) via a generating function namely 


P = 



1 \ 


G{xi) 


Xi=0 


(17) 


We hrst consider the case where all CKjS are positive. We hnd that the generating function 
G{xi) is given by 


G{x,) = 


where 




ni=i(i - ‘^biXi + (hf 




do — bo^ 

di = -bido + rCiCo = -bido - rcifo, 
dij = -bjdi + rCj{bifo + cdo) = -bjdi - rcjfi, 
dijk bj^d-ij TCk(bjfi “ 1 “ Cjd{^ = b^dij 

dijkl bidijk '^Cl(bkfij “ 1 “ ^kdij) = bldijk f'^lfijki 

dijklm bjYidijkl “t" "^CmiJjlfijk 


( 19 ) 


We see that the dij,,, coefficients have the structure of nested sums. 

So far we have considered the case a* > 0. When one or more of the OjS are negative 
we simply make the replacement in the generating function bj —)■ bj/{b‘^ — rc^) and Cj —)■ 
—Cj/{b‘^ — rc^), for every j with aj negative. 

We need to express (fT8|) in terms of kinematic data. This can be achieved without the 
need of knowing the explicit solutions of the scattering equations by considering special 
cases of ffT6D and then using our algorithm in order to evaluate these special cases. There 
are many ways to do this and here we provide one of them. We have 


&n = P 


(0,0,0,0,0)i 


bl — rcl = 


bi = 


-P(2,0,0,0,0)-P(-l,0,0,0,0) 


26o&i — P(i,o,o,o,o) 

-P(-1,0,0,0,0) 

^0 -P(1,0,0,0,0) 


2P(i,o,o,o,o)-P(-i,o,o,o,o) “ 26 q 


( 20 ) 
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and similarly for the remaining indices. For the nnmerator of fflSj) we have di = P(i,o,0,0,0) ~ 
2bobi, whereas in order to hnd du we can consider the 0,0,0) case, in order to hnd ^123 
the P(i,i,i,o,o) case and so on. Performing the algebra we hnd 


= - 1 -’ 

^2,2+1^2+2,2+3 



, 2+1 


1 

ki+2,i+i 




^i,i+2^i+1,1+4 
^i,i+l^*+2,i+4 


2bi 

-bibo + rCiCo 


^*+l,i+3^i+2,j+3 ^i,i+2 ^*+l,i+4 ^*,*+3^i+3,i+4 ^*,*+2 

^i,i+l^fi+2,i+4 ^i+2,j+4 ^i,i+l ^fi,'ii+l^i+2,i+4 ^*,*+1 

^i,i+2^i+l,i+4(^i,i+3^i+2,i+3 “b ^i,i+3^i+l,i+2 “b ^i+l,i+2^i+l,i+3) 
fcl2^23^34/i;45^15^i+2,i+4 


^i+l,i+4 

^i+2,i+4 


( 21 ) 

This is enongh data to determine the generating fnnction for all of the 2^ = 32 possibilities 
of signs of the different OfjS. From fllbp we see that we can simnltaneously change the signs 
of Co and CjS withont altering the hnal resnlt. This can also be seen from the fact that the 
coefficients of the generating fnnction always involve products of cq and CjS that cancel 
the sign. In fl?I|) we have given two equivalent expressions for the 6 jS. One can go from 
one to the other using conservation of momentum. One might think that since in flT^ 
we have nested sums, the coefficients will become increasingly complicated. This is not 
necessarily the case, since cancellations can and do occur. For example, in the case of all 
Qfj > 0 we simply have dijkim = 0. We should hnally mention that for the n = 5 case one 
can alternatively use the explicit solutions of the scattering equations in order to hnd the 
coefficients in fl2T|) . 

For completeness, let us briehy consider the n = 4 case. We have that the amplitudes 
are linear combinations of the following fundamental quantity 


roots 


1 

det'd* 


TT 1 f 0'13C’‘24A ^ / <713*^24^ 


( 22 ) 


Since we only have one root that we need to sum over, we can explicitly evaluate fl 22 |) to 
be 


- ^ + ^1 + 7^ 


12 


h 


23 


'423 


k 


12 


ai 


i + r^ 

^23 


a2 


(23) 


One can easily hnd a generating function that reproduces the result. For example, for 
cij > 0 we get 




( 24 ) 


We see that for Xi = 0 the generating function gives the expected answer for the color 
ordered 0 ^ theory. 


5 Examples 

As we have mentioned the amplitude is not always manifestly expressed as a linear com¬ 
bination of the fundamental quantity ([13]). Nevertheless, we can always bring it to the 
desired form and we demonstrate how this can be done in the case of Yang-Mills. 
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We choose to remove the hrst two rows and columns of the reduced pfaffian in (|1]). 
Then we are left with an 8 x 8 antisymmetric matrix, whose naive pfaffian expansion 
contains 105 terms. Most of the terms are already of the form f[T5]l . except for three cases. 
The first case involves one diagonal element of the matrix Cab, the second case involves a 
product of two diagonal elements of Cab, whereas the last case consists of three diagonal 
elements. In all of the three cases the treatment is the same. We pick up one of the four 
terms of each Caa and we replace it using conservation of momentum. Let us illustrate 
this with an example. Upon expanding the reduced pfaffian of the Yang-Mills, one of the 
105 terms has the form 


1 

C’'240'l5C’'45 


f ^3'kl ^ ^3‘^4 f3'^5 

V <^31 <732 <734 <735 


(25) 


We have omitted an overall factor that depends purely on helicities and momenta and 
we have kept only the part that contains the cr* variables. We now replace in the above 
expression the momentum ki —k 2 — to get 


1 

<724<715<745 


<7i2 

<7i3<723 


e3'^2 + 


<7i4 , 

- Gs'ki + 

<7i3<734 


<7i5 

<7i3<735 



(26) 


We have now achieved our goal. Combining all elements together the contribution of our 
term to the scattering amplitude is proportional to 


^ rlpr'd> n 


1 f <7i2<734 


roots 


det'<I) f \ crU+i V ^13^^24 


^3‘^2 ~ 


<7i4<723 

<7i3<724 


erki - 


<715<723<734 

<713<724<735 


^3'h 


(27) 


In terms of our fundamental quantity (|T^ . the above expression becomes 

e3'^2.P(-l,0,-1,0,0) — e3-fc4P(0,-l,0,0,0) “ ^3-k5P{0-1,-1,0-1)- 

The three Ps appearing in (1^ can be easily computed from (ITTjl . We find 

111 

+ + 


_ b2do - rcsCo 
t"(o,-i,o,o,o) - —^2 - 


ret 


^12^45 


k,ok: 


12^34 


^15^34 ’ 


(28) 


(29) 


where the denominator h\ — rc\ comes from the fact that the 0:2 is negative in our example. 
The same way we find 


-P(-i,o,-i,o,o) — 


and finally. 


(5i& 3 -F rciC3)do - (5iC3 -F hciYcQ 
{h\-rcl){hl-rcl) 


-P(0,-1,-1,0,-1) — 




1 1 
+ 


15«^23 


^23^45 


(30) 


(31) 


^ 12^45 

The case of the full Yang-Mills and gravity can be treated the same way, whereas for 
more general theories we expect similar considerations. 
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6 Discussion 


The main motivation of this work was to nicely organize and even calculate the tree am¬ 
plitudes of theories whose S-matrix can be described through the scattering equations. In 
doing so, we considered the most general quantity consistent with SL{2, C) invariance, 
that in general depends on conformal cross ratios of the variables that appear in the 
scattering equation. Then, all amplitudes can be written as linear combinations of this 
quantity, with coefficients that depend on kinematic data. We have found that our fun¬ 
damental quantity can be nicely expressed through a generating function, that we have 
explicitly calculated for the hrst non-trivial case, namely n = 5. Although the solutions 
of the scattering equations are complicated in nature, we have argued that knowledge 
of them is not necessary to evaluate the amplitude. We have also presented a simple 
argument why the amplitude is a rational function of the kinematic invariants. 

Although all hve point amplitudes can also be obtained by brute force, it might be time 
consuming to simplify all square roots appearing, whereas the hnal answer can be given 
in a disorganized form. Experience shows that for practical purposes the computation 
can be simplihed using the polynomial form of the scattering equations, but even that 
does not solve the problem of organization. Our proposed generating function gives the 
answer in a neat way. One can also see that it is rational in the kinematic invariants, it 
has the structure of nested sums and it knows about all signs coming from combinatorics. 

Our method can certainly be applied to the n = 6 case using our algorithm and 
the results of [^, whereas a generalization to the arbitrary n case remains to be seen. 
The general case is expected to be captured by a generating function of n{n — 3)/2 
variables, which is simply the number of all possible aijS with i < j after we subtract 
the n conditions coming from the SL{2, C) invariance of the problem. The form of our 
fundamental quantity flT^ can be easily extended to the general case and it involves the 
product of n{n — 3)/2 cross ratios each one appearing an integer number of times. The 
form of the generating function in the general case is also easy to be found. It will involve 
a denominator of n{n — 3)/2 polynomials, each one of degree (n — 3)!. The difficult 
part of the computation is to express the various coefficients of the generating function 
in terms of kinematic data. One possible way to achieve this would be to assume the 
form of the generating function and £x its coefficients by studying simple cases, where 
we already know the answer. It is also worth to investigate whether there is an even 
simpler but equivalent form of our generating function for n = 5, having always in our 
mind application to the general case. 
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